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Using the Relativistic Mean Field Theory (RMF) it is shown that different 
Q \ proton fraction which is directly connected with the neutron excess and 

with the asymmetry of the system affects proto-neutron stars parameters 
and changes their composition. The obtained form of the equation of state 
allows to construct the mass-radius relations and shows that the increasing 
J> ' asymmetry creates more compact stars. The inclusion of 5 meson together 

with nonlinear vector meson interaction terms and magnetic field make 
q ! this effect even stronger. 
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Introduction 



O 

Properties of dense matter in strong magnetic field has been the subject of in- 
vestigations in astrophysics of white dwarfs, proto-neutron and neutron stars. 
Observations of pulsars suggest large surface fields of the order of 10 14 G Jlfl. 
This very high value indicates the existence of even stronger interior magnetic 
fields. In fact the virial theorem suggests that interior magnetic fields can reach 
the value of the order of 10 18 G. The influence of this extraordinary high mag- 
netic field on neutron star matter properties are of particular importance after 
the discovery of magnetars [I]]. Their observations suggest that in the case of 
these objects we are dealing with young neutron stars with extremely strong 
surface magnetic fields ~ 10 15 G which in turn give an interior field of the 
order of 10 18 G. Also properties of proto-neutron stars, under the influence of 
magnetic field should be examined. Assumptions concerning composition and 
the equation of state of hot, lepton rich matter have been studied by many 
authors (Strobel et al. 1999, Bombaci et al. 1995, Takatsuka 1995, Mahka et 
al. 2001) 0J^J^,|]. Matter inside a proto-neutron star is highly degenerate and 
chemical potentials of its constituents are a few hundreds of MeV. Like in a 
neutron star the strength of magnetic field of a proto-neutron star changes 
from 10 8 G at the surface up to 10 18 G in the center ||. A proto-neutron star 
|7|J|U is a result of a supernova explosion which forms low entropy core with 
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trapped neutrinos. The core is surrounded by a low density, high entropy man- 
tle. The construction of a proto-neutron star model is based on various realistic 
equations of state and gives a general picture of proto-neutron star interiors. 
The more complete and realistic description of a proto-neutron star requires 
taking into consideration effects of finite temperature and nonzero magnetic 
fields. Using the relativistic mean-field theory approach the adequate form of 
the equation of state enlarged by contributions coming from magnetic field 
and temperature is constructed and serves as an input to the Oppenheimer- 
Volkoff-Tolman equations. The theory considered here comprises electrons, 
neutrinos, scalar, vector-scalar and vector-isovector mesons. The RMF theory 
implies that the nucleon interactions appear through the exchange of meson 



fields [y^,yj,|^,|l3|] . This theory is very useful in describing properties of nu- 
clear matter and finite nuclei. Its extrapolation to large charge asymmetry is 
of considerable interest in nuclear astrophysics and particulary in constructing 
proto-neutron and neutron star models where extreme conditions of isospin are 
realized. The model considered describes high isospin asymmetric matter and 
it has to be extended by the inclusion of isovector-scalar meson a (980) (the 
8 meson) and nonlinear vector meson interaction terms. Knowing the form of 
the equation of state is the decisive factor in determining properties of proto- 
neutron stars such as: central density, mass-radius relation, crust extent or the 
moment of inertia. The essential goal of this paper is to obtain the equation of 
state for proto-neutron star matter within the described above model in the 
relativistic mean field approach (RMF) and to study the influence of proton 
fraction on proto-neutron stars parameters. The variable proton fraction to- 
gether with the inclusion of 5 meson and nonlinear vector meson interactions 
alter the proto-neutron stars chemical composition. This in turn affects the 
properties of the star making the proto-neutron star more compact. 



The generalized Relativistic Mean Field approach 



The Lagrangian function for the system can be written as a sum of a bary- 
onic part Cb which includes baryon-meson interaction terms, mesonic part 
Cm containing additional interactions between mesons which mathematically 
express themselves as supplementary, nonlinear terms in the Lagrangian func- 
tion, the leptonic part Cl, the Lagrangian density function of the QED theory 
£qed and the gravitational term Lq 

£ = f~-B + C-L + C-M + C-G + £-QED- (l) 

The final form of the Lagrangian function is given by 
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C = iip-fD^ - tfj{M N - g a cp - I 3N g s r a 5 a )tfj - K^a^F^ (2) 

+i £ LffD^Lf - £ 9f(L f He Rf + h.c) - \d^d»y - U{ip) 
/=i /=i 



with the covariant derivatives -D^ and defined as 



D^ = d^ + -ig p p°c a + igjujp + iQA^, D^ = d^ + iQA^. (3) 

In the case of nucleons Q takes the value e for protons (e is the electron charge) 
and for neutrons (Q = e, 0). In the leptonic sector (Q = — e, 0) where — e is 
given for electrons and muons and for neutrinos. H is Higgs field and this 
field has only the residual form 




where the value v = 250 MeV comes from the electroweak interaction scale. 
In a strong magnetic field, contributions coming from the anomalous magnetic 
moments of protons and neutrons (k^, N = {p, n}) have to be also consid- 
ered. The anomalous magnetic moments introduced via the minimal coupling 
of nucleons to the electromagnetic field tensor can be represented phenomeno- 
logically by interactions of the type Kn^ct^F^i/j , where a pv = §{7 M , 7^} fyfl . 
Considering both the physical conditions and chemical equilibrium which are 
likely to be satisfied in a proto-neutron star one can construct a proto-neutron 
star model with a temperature T equals 20 MeV and magnetic field strength 
B ~ 10 3 B^ where B e c is the value of the critical magnetic field for an electron 
B e c = m 2 e j | e |= 4.414 x 10 13 G. The model examined in this paper represents 
the one with moderate value of magnetic field. In this case the contributions 
to the equation of state coming from anomalous magnetic moments are small 
(~ lMeV) in comparison with the energy of magnetic field. Thus this rela- 
tively small value of magnetic field allows to neglect the anomalous magnetic 
moments unlike to the case describe in the paper by Prakash et al. [fOJ where 
authors consider the case of extremely strong magnetic field of the order of 
10 18 G. Realistic proto-neutron star models describe electrically neutral, hot, 
high density matter being in (3 equilibrium. The last condition implies the 
presence of leptons which is expressed by adding the Lagrangian of leptons to 
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Table la 

The parameter sets of the model [p^j|To| . 





9 P 


98 


9a 


9u 


M a (MeV) 


TMl 


9.264 





10.0289 


12.6139 


511.12 


TMl + nonl. 


10.875 


3.5 


10.0289 


12.6139 


511.12 


GM3 


8.5417 





7.1857 


8.7041 


450.00 



Table lb 

The self-interacting coupling constants |14] , |15 





92 (MeV) 


#3 


c 3 


A v 


A 4 


c 


TMl 


1427.18 


0.6183 


71.3075 








TMl + nonl. 


1427.18 


0.6183 


71.3075 


0.008 


0.001 


0.5 


GM3 


3016.52 


-6.4546 












the Lagrangian function. The fermion fields which are included in this model 
consists of neutrons, protons, electrons, muons and neutrinos 
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As a proto-neutron star matter is of sizeable asymmetry the additional 5 
meson has been included and the meson sector is composed of isoscalar (scalar 
a and vector uj) and isovector (scalar S and vector p) mesons. The potential 
function U((p) = ^M^ip 2 +^2^ +\gz^ has a very well known form introduced 
by Boguta and Bodmer ||16|| . Nucleon masses are denoted by M N (N=p,n) 



whereas M w , M p , M a and Ms are masses assigned to the meson fields. They 
are taken at their experimentally values: M^ = 783 Me V, M p = 770 MeV 
and Mg = 980 MeV. The nonlinear vector meson interaction terms have been 
added in order to give a detailed account of a matter inside a proto-neutron 
star. This has been done through the construction of adequate form of the 
equation of state. To determine the equation of state at finite temperature and 
nonvanishing magnetic field some parameterizations have been used. The first 
one is known as TMl parameter set |15|, the second one noted as TMl + nonl. 
includes nonlinear terms in the meson sector. The third is the GM3 parameter 
set [|14|]. The parameters entering the Lagrangian function are the coupling 
constants g p , gs, g a , gw for meson fields and self-interacting coupling constants 
52 #3, C3, Ay, A 4 and (. All these parameters have been chosen to reproduce 
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Table 2 

The nuclear matter properties. 





E b (MeV) 


m eff /M 


n s (fm~ 3 ) 


K (MeV) 


a sym (MeV) 


TM1 


-16.3 


0.658 


0.1455 


281.99 


36.82 


TM1 + nonl. 


-16.3 


0.664 


0.1455 


281.99 


32.22 


GM3 


-16.3 


0.78 


0.153 


240 


32.5 



properties of symmetric nuclear matter at saturation and they are collected in 
Tables la and lb. In the original GM3 and TM\ approaches neutrinos are not 
included since these parameterizations describe a neutron star matter. In this 
case neutrinos leave the star unhindered. In order to construct the neutron 
star model through the entire density span the addition of the equations of 
state, characteristic for the inner and outer core, relevant for lower densities 
is necessary. In the result the composite equation of state for which the TM1 
or GM3 parameter groups describing the neutron star core are supplemented 
with the Negele-Voutherin (NV) + Bonn describing inner crust can be used. 
The equilibrium properties of nuclear matter which are known with reasonable 
precision are: the nuclear saturation density ns, the binding energy E^, the 
incompressibility K , the symmetry energy a sym and the effective nucleon 
mass. Nuclear matter properties calculated for given parameter sets (Tables 
la and lb) are summarized in Table 2. The presence of additional nonlinear 
p meson interaction terms and 5 meson which carries izospin contribute to 
the symmetry energy. Thus the value of the parameter g p has to be redefined 
for the TM1 + nonl. parameterization. Parameters g p , gs, A v , A 4 and ( are 
not chosen arbitrary they have to be given in a specified form adequate to 
reproduce the value of the symmetry energy a sym . The vector potential for 
the electromagnetic stress tensor F pv is given by A p = {A = 0, A^} where 

Ai = ~£ilm% . 

The symmetry in which uniform magnetic field Bo lies along the z-axis has 
been chosen chosen B™ = (0,0, B z ). Considering the influence of an external 
magnetic field on proto-neutron star constituents one can start with determin- 
ing the dispersion relations for them. The dispersion relation for an electron 
in a magnetic field is 



e e = \Jp 2 z + m 2 e + 2neB z , 
where n is the Landau level, p z is the electron momentum along the z-axis and 



m e is the rest mass of the electron 17 . Along the field direction a particle 



motion is free and quasi-one-dimensional with the modified density of states 
which is given by a sum 
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5 n0 denotes the Kronecker delta [17]]. The spin degeneracy equals 1 for the 
ground (n = 0) Landau level and 2 for n > 1. Similar results can be obtained 
for muons by replacing the electron quantities by the corresponding muon 
quantities. The energy dispersion relation for protons for arbitrary Landau 
level in a magnetic field is JT^] 

£ P ,s = \jp 2 z + (y/M^ff + 2nQB z + sk p B z ) 2 (5) 

and s = ±1 being "spin" number. For moderate values of magnetic fields 
strength this equation reduces to the following one 



e p ~y/p> + Aq teff + 2neB z . 

A neutron can interact with an external electromagnetic field by means of the 
Pauli non-minimal coupling. Then the energy dispersion relation for neutrons 
in a magnetic field is given by 

Sn,s = ^p\ + ^Ml eff ^p{+SK n B z f (6) 

where m = p z and p± are the components of the neutron momentum parallel 
and perpendicular to the magnetic field. As it is usually assumed in quantum 
hadrodynamics the mean field approximation is adopted and for the ground 
state of homogeneous infinite matter quantum fields operators are replaced 
by their classical expectation values. One can separate mesonic fields into 
classical mean field values and quantum fluctuations which are not included 
in the ground state and thus mesons fields are replaced by their mean values 
a =< >, d a =< 5 a >= 5 a ' 3 d , Wfl =< >= S^qWq and r£ =< p£ >= 
S^S^o r . The Dirac equation at the mean field level for the nucleon quasi- 
particle has the form 

(i 7 ^ - M N>eff - K N a^F^)ij = (7) 

with Mjv, e // being the effective nucleon mass generated by the nucleon and 
scalar fields interactions and defined as 



M Pteff = M5 P = M- g s a - g s d , (8) 
M n>eff = M6 n = M - g s a + g s d . (9) 

The main effect of the inclusion of 5 meson becomes evident studying prop- 
erties of proto-neutron star matter especially nucleon mass splitting and the 
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form of the equation of state. With the effective nucleon mass Mjv, e // one can 
redefine the proton and neutron chemical potentials 



/j p = e p + g^wo + -g p r (10) 
Hn = £n + 9u>w Q - -g p r (11) 

where e p and e n is given by relations (§) and (|5|) . The obtained effective nucleon 
chemical potentials allow to define conditions that are regarded as essential 
to uniquely determine the complete equilibrium composition of proto-neutron 
star matter at given density. These conditions arise from charge neutrality, 
baryon and lepton number conservation. The last one is strictly connected 
with the assumption that proto-neutron star matter is opaque to neutrinos. 
Neutrinos are trapped inside the matter and this through the requirement of (3 
equilibrium puts certain limits on the chemical composition of a proto-neutron 
star matter. For the case to be considered the processes of /3-decay and inverse 
/3-decay take place |18| 

p + e <-> n + v e . (12) 

For large enough electron Fermi energy it is energetically favorable for elec- 
trons to convert to muons 

fi + v e <-> e + z/ M . (13) 

The chemical equilibrium established by the above processes is given by re- 
lations between chemical potentials of proto-neutron star constituents and 
impose constrains on them. 



/V = yU e + V P ~ (14) 
= /V -He + Hp- (15) 



From the equations (|10D and ([11]) the electron neutrino chemical potential can 
be expressed in the following way 



f*v e = He + £p - £n + 9pT0- (16) 

With the assumption that only electron neutrinos are captured in the star 
core (//j. = 0) the following relations for electron neutrino and muon chemical 
potentials can be obtained 



/V = -fJ>v e + He = -Ep + £n~ 9 p r - (17) 



7 



Both electron neutrino and muon chemical potentials depend on nucleon asym- 
metry. The muon chemical potential depends on the asymmetry of the whole 
system through the value of the p meson field which is directly connected 
with the neutron excess and the difference of nucleon energies. Comparison 
between the matter with high asymmetry (Y p = 0.11) and the one with low 
asymmetry indicates that in the former case the muon chemical potential is 
greater than that of the later one. This translates to the lower number of 
muons in the proto-neutron star matter with the proton number Y p of the 
order of 0.38. The increase of neutron excess is connected with the increasing 
number of muons. Their contribution to the pressure and energy density of the 
system become significant. To construct the equation of state of the system 
the energy-momentum tensor has to be calculated. The pressure is related to 
the statistical average of the trace of the spatial components of the energy- 
momentum tensor P = | < Ta > whereas the energy density e equals < T 00 >. 
In general in the presence of magnetic field the pressure can be written as a 
sum of the isotropic part which includes contributions coming from electrons, 
nucleons, mesons and the electromagnetic anisotropic part T^ u is given by 



In cartesian coordinates the electromagnetic part of the energy momentum 
tensor in the flat space-time has the form 



I ' \B 2 
\B 2 



■ in; 








\ 



\B 2 







-\B 2 



(18) 



On the other hand in polar coordinates T^ v has been written as follows 



' \B 2 N 

-\B 2 cos26 \B 2 rsm26 

±B 2 r sin 20 \B 2 r 2 cos26 
v \B 2 r 2 sm 2 9 ) 

Taking the average over all directions allows to obtain the isotropic form of the 
energy momentum tensor relevant for the spherical symmetric configuration 
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1 \B 2 
i5 2 

6 

i5 2 

V o o o \b 2 } 

The assumption of the positivity of the total pressure P has been made. Any 
negative contributions to the pressure diminish it and lead to the decrease 
of a star radius. The described above averaging approximates the star as an 
spherically symmetric object with the radius R which is smaller than the 
equatorial radius of an anisotropic star. (The case of an anisotropic star is 
calculated in details in the paper by Konno et al. The condition of 

the positivity of the total pressure establishes a limiting value of the density 
noted as p cr u and at the same time puts limit on the spherical symmetric 
approximation (the deformation of the star is neglected) . For densities smaller 
than the critical one the total pressure became negative and there is no gravity 
compensation in this direction. The total pressure and the energy density 
of the system supplemented by corrections coming from magnetic field are 
obtained by adding contributions from mesons and fermions 



P = P f + Pqed ~ U(a) + ^M 2 r 2 - ±M 2 w 2 - l -M 2 d\ + A^r> 2 + 
+<7^A 4 r V + ^c 3 w 4 + ^( 9 y 



e = e f + sqed + U(a) + -c 3 w 4 + -M> 2 + -M 2 r 2 + 3A v g 2 g 2 p r 2 w 2 + g 2 a g 2 p 



3 

T" ° ' 2 ~ " ' 2 
A 4 r V 2 + \C9 A A + \M 2 A 



where / = (i,n), i = (e,//,p) and e QED = \B 2 ,P QED = \B 2 . The fermion 
pressure is defined as 



T 77? 4 °° 

p i = TT E EI 2 - W-iA+CVti, S? + 27en) 

47r s n =0 



where the Fermi integral 



f (a + x 2 ) A/2 x r? da; . . 

h, v ±(u,a) = / 1 , 7 =4- , 19 

was used, ^ = /lo/m, U = k B T /mi, m = Zi/t h 7 e = B z /B e c and q = 
m^eff/rrii. The last quantity equals 1 for electrons and muons. The critical 
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Figure 1. The lepton chemical potentials (e, fi, v e ) as a function of the baryon density. 

magnetic field strength for electron equals Bl=m 2 e / \ e |= 4.414 x 10 13 G. The 
fermion energy density is also defined with the use of the Fermi integral 



47r s n=0 

+Ii fi ,-(zi/ti,<$ + 2i e n)}. 



Conclusion 



Analyzing properties of proto-neutron star matter some features which are 
relevant for proto-neutron stars have to be considered. It is a necessary as- 
sumption that the matter is opaque to neutrinos. As a result neutrinos are 
trapped dynamically inside the matter on the diffusion time scale and in con- 
tradiction to neutron star matter the total lepton number Yl = Yg+Y^+Y^ is 
kept constant. Theoretical predictions determine the value of Yl at the onset 
of trapping at about 0.35-0.4. Since a proto-neutron star is formed during a su- 
pernova explosion it is safely to make an assumption about finite temperature 
in the interior of a star. This is limited by the requirement of constant value of 
entropy in a proto-neutron star core. In this paper the temperature is fixed and 
equals 20 MeV. Thus conditions which have to be considered namely: finite 
temperature, lepton number conservation, (3- equilibrium, charge neutrality 
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Figure 2. The equation of state for the proto-neutron star for different values of 
the proton concentration Y p and strength of magnetic field. TM1 model with the 
nonlinear terms. 




R (km) 



Figure 3. The mass radius relations for different proton fractions Y p (T = 20, MeV ', 
B ~ 10 3 B^., 5^=4.414 x 10 13 ). For comparison the mass radius relation for the 
proto-neutron star without magnetic field (B = 0) is also presented). 
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Figure 4. The mass radius relations for neutron and proto-neutron stars for different 
parameter sets. 

and additionally nonzero magnetic field allow to construct a realistic proto- 
neutron star model. All of these mentioned above conditions influence medium 
properties such as nucleon effective masses, density and proton-neutron asym- 
metry and in consequence constrain the nucleon and lepton species inside the 
star core. The proton fraction Y p = n p J (n p + n n ) which is determined by (3 
equilibrium, for matter transparent to neutrinos (neutron star matter) takes 
the value of the order of ~ 0.1 whereas when neutrinos do contribute they 
displace the equilibrium proton fraction to the higher value Y p ~ 0.3 — 0.36. 
For the purpose of comparison all calculations have been performed for de- 
creasing value of proton fractions Y p = 0.34,1^, = 0.22, Yp = 0.11, respec- 
tively. This changing proton-neutron asymmetry through relation (|T7|) alters 
the lepton chemical potentials. Figure 1 displays chemical potentials of elec- 
tron neutrinos, electrons and muons for two extreme values of Y p (Y p = 0.11 
and Y p = 0.34). They have been examined for moderate strength of magnetic 
field. If the proton fraction equals Y p = 0.34 (the proto-neutron star matter) 
electron and neutrino electron chemical potentials are increasing functions of 
a density. Neutrinos are trapped inside the matter and their chemical poten- 
tial differs from zero (see relation |l6|). For the proton fraction Y p = 0.11 the 
matter inside the star resembles that of a neutron star which is characterized 
by larger value of asymmetry. Neutrinos are about to leave the system and 
their chemical potential vanishes \i v — > 0. This picture shows the profile of 
neutrino chemical potential in the core of the star for the chosen values of 
Y p and for chosen parameter set (TMl+nonl.). For the star with the value of 
asymmetry typical for a proto-neutron star matter (Y p ~ 0.38) muons are of 
insignificant meaning. However, they become much more important for lower 
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value of proton number. The asymmetry entered into calculations at different 
levels starting from the displacement of proton and neutron chemical poten- 
tials fi p 7^ fi n . In a model with 5 meson exchange additionally M p ^ M n . There 
are also contributions coming from nonlinear vector meson interactions. The 
nuclear asymmetry diminishes the neutrino chemical potential and thus low- 
ers neutrino density. The influence of the asymmetry is most distinctive in the 
high density region inside a proto-neutron star. The forms of neutrino profiles 
in the proto-neutron star core depend on the asymmetry of the whole system. 
The behavior of chemical potentials translates to the forms of the equation of 
state and successively to the mass-radius relations. Results are obtained for 
two particular cases namely for the nonmagnetic one (B = 0) and for the value 
of magnetic field B ~ 10 3 5^. The presence of magnetic field makes the equa- 
tion of state stiffer (Fig. 2), whereas for more asymmetric matter the obtained 
equation of state is softer. Dots represent the limit of making the approxi- 
mation of spherically symmetric star in the presence of magnetic field. For 
more asymmetric matter the dot is shifted towards higher densities. In Fig. 3 
one can compare the influence of changing proton fraction on star parameters. 
The obtained sequences of mass-radius for nonmagnetic and magnetic proto- 
neutron stars show that the increasing asymmetry diminishes proto-neutron 
star radii causing that the star parameters start to resemble that of neutron 
stars. This effect is even stronger in the presence of magnetic field. There are 
relatively minor changes in the value of maximal masses caused by the vari- 
ation of proton fraction. However, there is a tendency of the increase of the 
maximal mass for growing asymmetry. The star parameters are also sensitive 
to the presence of nonlinear vector meson interaction terms in the Lagrangian 
function (|j). This is clearly visible in Fig. 4. Comparison of results obtained 
for original GM3 and TM1 equations of state (for neutron stars) with that ob- 
tained for TMl+nonl. shows that the last one gives the configuration with the 
bigger value of maximal mass. In the range of lower densities the TMl+nonl. 
parameter set gives the star with the smaller radius. Thus in the case of high 
asymmetry nonlinearities create more compact objects. Nonlinearities alter 
proto-neutron star parameters insignificantly. This is caused by lower value of 
the asymmetry. 
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